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Abstract 



When two perfectly conducting inclusions are located closely to each other, the electric 
field concentrates in a narrow region in between two inclusions, and becomes arbitrarily large 
as the distance between two inclusions tends to zero. The purpose of this paper is to derive an 
asymptotic formula of the concentration which completely characterizes the singular behavior 
of the electric field, when inclusions are balls of the same radii in three dimensions. 
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1 Introduction and statement of results 

Let Di and Dy, be bounded, simply connected and convex domains in Mr, d = 2,3. Suppose that 
the conductivity of the inclusions is oo, in other words, inclusions are perfect conductors. We 
consider the following conductivity problem: 

Au = in R d \ £>i U£> 2 , 

u = Cj ; (constant) on dDj, j = 1, 2, (1.1) 

u(x) - H(x) = 0{\x\ 1 - d ) as|x|->oo, 



where H is a given harmonic function in R d so that — VH is the background electric held in the 
absence of the inclusions. The constant value C,- on dDj is determined by the condition 

in . 

O ■ r du 

C*V / T777T dCT = ° fori = 1,2. (1.2) 



X 



JdD] dvti) 

Here and throughout this paper u^> is the outward unit normal to dDj. 

The gradient of the solution Vw represents the electric field (with the opposite sign) in the pres- 
ence of inclusions and the stress field in two dimensional anti-plane elasticity, and it may become 
arbitrarily large as the distance between two inclusions tends to 0. It has been proved that the 
generic rate of the gradient blow-up is e -1 / 2 in two dimensions [3J |H [5J |H1 EJ HH [H] and leloge^ 1 
in three dimensions [SJ 1121 113) . where e is the distance between two inclusions. Occurrence of 
the gradient blow-up depends on the background potential (the harmonic function H in (|l.l|l ) and 
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those background potentials which actually make the gradient blow up are characterized in [3] 
when D\ and D2 are disks. 

The results mentioned above are estimates of the gradient of the solution from above and below, 
namely, 

< |V«| < -^r + C 3 (1.3) 

for some positive constants C\ . C2 and C3 where 

if d=2, 

m = { ;,„„i ifd=3 . d-4) 

The constants C\ and C2 can possibly be depending on the background potential H. 

The interest of this paper lies in the asymptotic behavior of V« as the distance between two 
inclusions tends to 0. Since the singular behavior of Vu occurs in the narrow region in between 
two inclusions, we arc particularly interested in its behavior there. In this regards, a complete 
characterization of the singular behavior of Vu has been obtained when inclusions arc disks |10| 
and strictly convex domains in M 2 [T] . Let D\ and D2 be disks in R 2 of radii t\ and , respectively, 
and let Rj be the reflection with respect to dDj, j = 1,2. Then the combined reflections R1R2 
and R2R1 have unique fixed points, say fi £ D\ and f2 G D2. Let 

/ l (x) = i-(log|x-f 1 |-log|x-f 2 |) (1.5) 

(see section [2] for a discussion on the function h). It has been proved that the solution u to (|1.1D 
can be expressed as 

ri + r 2 

where c is the middle point of the shortest line segment connecting dD\ and dD2, n is the unit 
vector in the direction of fa — fi, and |Vp(x)| is bounded independently of e on any bounded subset 
of R 2 \ (D\ U L>2). So the singular behavior of Vu is completely characterized by V/i. In particular, 
it can be shown using (jl.6p that the maximal concentration of Vu occurs along the shortest line 
segment connecting dD\ and dD2, and on that segment 

Vti«^,&n.Vff)(c)a (1.7) 

A complete characterization of the gradient blow-up like (|1.6[) has been obtained in [T] in the 
case when inclusions are strictly convex domains in R 2 by using disks osculating to convex domains. 
It is worth mentioning that the stress concentration factor for the p-Laplacian was derived in [5] . 

The purpose of this paper is to derive an asymptotic formula of Vu which characterizes its 
singular behavior when D\ and D% are balls of the same radii in three dimensions. 

In order to state the main result of this paper in a precise manner, let us fix notation. Let D\ 
and D2 be balls of radius r in three dimensions and ci and C2 their centers. Let c be the middle 
point of Ci and c 2 , and n the unit vector in the direction of C2 — Ci, i.e., 

Ci + C 2 C2 — Ci 



2 ' |ca-ci| 
Let Rj, j = 1, 2, be the reflection with respect to dDj, i.e., 

„ , . r(x — cA 

R^) = j^+cj, 
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and let, for k = 0, 1, . . ., 

ip 2 k = (R 2 Rl) k C 2 , (1 

\ P2fe+ i = R 2 {RiR 2 ) k c 1 . 

We emphasize that p„ is contained in D 2 and monotonically converges to p as n —> oo where p is 
the fixed point of the combined reflection R 2 R\. Let 

Mn=, " r, n=l,2,..., (1.9) 

|ci - p„| 

and 

n 

q a = 1 and q n = ]^[ ^, n > 1. (1-10) 

Let p(x) be the distance from x to the line connecting Ci and c 2 , i.e., 

p(x) = |(x-c)-(x-c,n)n|. (1.11) 

The following is the main result of this paper. 

Theorem 1.1 Suppose that the radius of the balls is much larger than the distance between them, 
i.e., e <C r. The gradient Vu of the solution to (II. ip can be expressed as 

Vu ( x )=n M^l ? . 2 . (n + r ? (x))+V g (x) if p( x ) < (1.12) 

|loge| (e + rp(x) ) [log e| 

where 

oo 

C$j = 2J2 ^ (H( Pn ) - H(-p n )) , (1.13) 

n=0 

|V<7| is bounded on any bounded region in R 3 \ (Z?i U D 2 ) regardless of e, and 

|r?(x)| ^Clloger 1 (1.14) 

for some constant C > independent of e. 

Some remarks on Theorem 11.11 are in order. We first observe that the set p(x) < r|loge|~ 2 
where (|1.12JI holds is a narrow region in between two spheres. The formula (|1.12[) shows that the 
major singular term of Vit is in the direction of n, and that if p(x) = constant, then intensity of 
the field is constant. Note that the level set where p(x) is constant is a cylinder around the line 
connecting centers of two spheres. So the intensity of the field decreases radially from the line 
connecting two centers of spheres. The highest concentration of the field occurs when p(x) = 0, in 
other words, when x is on the line segment connecting two closest points on the spheres, and on 
the segment, 

Vw «-^n. (1.15) 
e|loge| 

Note that depends on e since p„ and q n do. The following theorem reveals the limiting behavior 
of C € H as e — >• 0. 

Theorem 1.2 We have 

C e H = C H + (Ve| log e|) as e -> (1.16) 

where 



' n 

n=0 



In particular, if /?(x) = 0, then 



lime|loge||Vu(x)| = \C H \- (1-18) 
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We emphasize that the occurrence of the gradient blow-up depends on the constant Ch'- if 
Ch 0, then it occurs. If Ch = 0, then either |Vu| is bounded or the blow-up rate is weaker than 
the generic rate (e| logeQ -1 . One can show for example that if the centers of the balls lie on the 
x-axis and their middle point is (0, 0, 0), and if H (x, y, z) = x 3 — 3xy 2 , then Ch ^ and hence | Vu| 
blows up as e — > 0. It is interesting to observe that this is in contrast with two dimensional circular 
case. In view of (|1.7[) . the blow-up occurs only when (n • ViJ)(0, 0)^0 (assuming c = (0, 0)). So, 
Vufi, y) blows up in two dimensions only when the background potential H has the linear term 
n • x. 

The main ingredient in deriving (|1.12[) is the singular function h which is the solution to 

Ah = in R d \ D 1 UD 2 , 

h = constant on dDj, j = 1,2, 

dh ,«. -(-!)*» , = 1,2, ( " 9 > 



{ h(x) = 0(|x| 1 " <i ) as |x| ->• oo. 

Such a solution exists and is unique (see [TJ E]). We emphasize that the constant values of h on 
dDi and on dD 2 are different, and because of that the gradient of h becomes arbitrarily large if 
the distance between D\ and D 2 is small. This function characterizes the singular behavior of the 
solution to (jl.l[) . In fact, if wc define the function g by 

u(x) = l\ dD2 ~l\ dDl h( X )+g( X ), xel^fAU^), (1.20) 

then one can see that g is harmonic in M d \ D\ U D 2 and g\dD ± = g\dD 2 i m other words, there is 
no potential difference of g on dD\ and dD 2 . So it can be shown in the same way as in [TU] that 
|V.g| is bounded on bounded subsets of M. d \ (Di U D 2 ). It means that the singular behavior of Vu 
is completely determined by ^° 2 _^J^° 1 V/i(x). Moreover, it is proved in [THrTS] that 

f dh f dh 

u\dD 1 -u\ aD2 = H— Tr -d(T+ H—p-da, (1.21) 
JdD 1 ouW J 9D2 <9/A 2 J 

which means that the potential difference of u is determined by the singular function h (and the 
background potential H). 

The function h was first introduced in |14j and used in a crucial way to derive estimates for 
the gradient blow-up in [T31 HH [T3]. It is worth mentioning that ( g ^ 1} , j^y) is an eigenvector 
corresponding to the eigenvalue 1/2 of the Neumann-Poincare operator associated with the interface 
problem (|1.1[) as shown in [TJ [7] . 

If Di and D 2 are disks, then h is given by (|1.5[) . In fact, dD\ and dD 2 are the Apollonian 
circles of the fixed points fi and f 2 , and hence |x — fi|/|x — f 2 1 is constant on dD\ and dD 2 . It 
is worth emphasizing that here the radii of disks may be different. If D\ and D 2 are spheres, 
it is proved in |13) that h is given by a weighted sum of the difference of the point charges: let 
r(x) = |^|x| -1 , the fundamental solution of the Laplacian in three dimensions. Then the singular 
function h is given by 



1 



/i(x) - =5^ — V q n (r(x - P „) - r(x + Pn )) . (1.22) 



n=0 1i 



=0 



This formula has been used in [T3| to derive estimates like (|1.3[) . We emphasize that in [T3] an 
upper bound for h is derived in a more general case when the radii of spheres are allowed to be 
different. In this paper we derive finer estimates of h for the purpose of deriving (|1.12p . 

This paper is organized as follows. In section [5J we review the construction of the singular 
function in [13] . In section [31 we prove some technical lemmas which are required to estimate the 
singular function. In section SI we derive an asymptotic formula of the singular function. In the 
last section, we prove Theorem 11.11 and Theorem 11.21 
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2 Singular functions on spheres 

Since the radius r is much larger than e, we may assume after scaling if necessary that r = 1. 
We may also assume the centers are on the a;-axis and c = (0, 0, 0) after rotation and shifting if 
necessary. We assume so in the sequel. It is also convenient to write e = 25 so that ci = (—1—6, 0, 0) 
and C2 = (1 + 6, 0, 0). Then, the function p defined in (jl.lip becomes 



p(x,y,z) = vV + ^ 2 , (2.1) 
and n = (1,0,0). Note that p„ defined by (|1.8[) satisfies 



P2k = (Ml) fc C 2 = -{R 1 R 2 ) k C 1 , 

P2fc+1 = — Rl{R-2Rl) k C2 = i? 2 (i?l-R 2 ) fc Ci. 



Define the function h\ by 



k=0 



hl{x) := V - , q2k+1 , • (2.3) 

^ \ |x + p 2fe | X-pafc+i / 



|x-p 2 fc+l| |ci - p 2 fc+l| |x - i?l(p 2 fe+l)| |x + p 2 fe+2| 

if x G dD\ , and 

<?2fc 52A: 1 <72fc+l 



Moreover, since 



we have 



Define h 2 by 



i /■ dhi, ^ i /• ^ 



fc=0 " ^ 2 A =0 



^ V |x-p 2fe | |x + p 2fc+ i| I 



(2.2) 



Then is harmonic in R J \ Dj U D 2 . Since the circle of Apollonius implies 

ly-Cj-Hx-ii^y)! = |x-y| for |y - c s \ > 1, x e 8D jt j = 1, 2, (2.4) 

we have 

<? 2 fc+l <?2fc+l 1 <?2fc+2 



|x + p 2fc | |c 2 + p 2fe | |x - i? 2 (-p 2fc )| |x-p 2 fc+l| 

if x G dD 2 . So we have 

hi\dD! = 1, hi\dD 2 = 0. (2.5) 



(2.6) 



Then h 2 is harmonic in ffi \ D\ U D 2 , and one can show similarly that 

h 2 \a Dl =0, h 2 \ dD2 = l, (2.8) 

and 

1 f 9h-) ^ 1 f dh 2 , _ N 

^Jo Dl dv ^ AttJ 9D2 dv ^ 

It then follows from (|23|) . (|2^|) . (|2T5]| . and (|23|) that the solution to (Ti~T9")) is given by 

1 1 °° / 1 1 \ 

fc(x):=-— =5 (/n(x) - /i 2 (x)) --^35 E<M I T~T ~ I — i ■ (2-10) 

47rE„=o9« V 7 47r}J n=0 g n ^ VI X ~P«I |x + p„|/ 
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Thus we have (|1.22|) . We also have 

h\dD 2 — h\g Dl 



(2.11) 



In the next section we derive fine properties of the sequences p„ and q n , which are used in 
deriving an asymptotic formula for h. 

3 Properties of the sequences p„ and q n 

Let p = (p, 0, 0) be the fixed point of the combined reflection R2R1 as before. Then one can easily 
see that p satisfies 

p = -— 4- \-l + 6, 



l + S+p 



so that 



p=V2S + 0(8) as £->0. 
Let p„ = (jp n , 0, 0). Then, pq = 1 + 6 and p n satisfies the recursive relations 

1 



One can further see that 



where 



Pn+l 



Pn = P 



l + S+p,-, 



1 + 6, n = 0,l,.. 



1 = P 



A 



n+l 



1 



A n+1 _ l J ^„ +1 _ j 

A . ±+s+p 



1 + 6-p 
Note that 

A = 1 + 2p + 0(6) = 1 + 2V25 + 0(6). 
In particular, the sequence p n is decreasing and converges to p as n — > 00 . 

1 1 



Since 



Cl - Pr, 



we have 



1 + 6 + p 
1 



= (1 + 6 -Pn+l), 

q n = (l + 6-p n+1 )q„ 



l + 5 + p n 

For a given 6 > 0, let N = N (6), N = N(6) and iVi = ATi(<5) be as follows: 



iV (5) = [|lo g 5|], iV(<5) = 



1 



V5. 



1 



tf|log<*| 



(3.1) 
(3.2) 

(3.3) 

(3.4) 
(3.5) 

(3.6) 
(3.7) 

(3.8) 



Here [•] is the Gaussian bracket. We use this notation for the rest of this paper. Since 6 is 
sufficiently small, we have 

N (6) < N(6) < iVi(tf). 
The following lemma was obtained in |13j . 
Lemma 3.1 There is a constant C independent of S such that 



n+l 



and 



forn < N(6). 



n + l 



C 



< cVs 



\Pn "Pn+l I < —n 



(3.9) 
(3.10) 
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We prove the following lemma. 
Lemma 3.2 Let N = N{5) and Nt = Nx(S) as before. 
(i) There is a positive C independent of S such that 

N 



y2 q n - Y] - 

n— n— 1 



< C and q n < C. 



(ii) p n -p> 2\flA~ n for all n. 

(iii) There is a constant C such that 



C 

Pn - P > — 



for all n< N. 
(iv) <p Nl -p < e -V(V*|iog*D. 
Proof. Since p n decays to p, we have from (13.71) 

q n < (1 + 6 - p) n ~ m q m for all n>m> 1. 
So, it follows from (|3.9[) that 



OO OO y \ oo 

n=AT n—N ^ ' n=N 



and 



A' 



< 



1 

E ?n - v - 

n— n— 1 

This proves (i). 

We have from (13.31) that for each n €E N, 



-N 



2p 



An+l -1 



So, (ii) follows from (|3.1[) . 

Now, suppose that n < N. Since A < 1 + 3p, using the inequality 

(1 + s) n < 1 + ns + ^n 2 s 2 (l + s) n 

which holds for all s > 0, we obtain 

A" < (1 + 3p) n < 1 + 3np + |nV(l + 3p) n . 

Since < TVp < 2 and (1 + i) 1 /' increases to e as i — > 0+, we have 

3np 



(l + 3p) n < l(l + 3p)*J " <e 6 , 

and hence, from the second inequality in p,14j) 

A n < 1 + C*np 



for some constant C independent of n < N and S. We then infer from f|3 . 3[) that 
Now, if n = Ni , then we have 



log(A n ) = nlog A > -i ^ > 



log<S|' 
and hence 

Now (iv) follows from (|3.13[) . This completes the proof. □ 
Lemma 13.21 (i) yields 

J2qn = ^\log5\ + 0(l). (3.15) 

n=0 

The following lemma provides the finer properties of p n and q n that are crucial in proving the 
main result of this paper. 



Lemma 3.3 (i) If N (S) <n< Ni(6), then 

q n l + Odlog*)- 1 : 



as S -> 0, (3.16) 



Pn - Pn+l y/pl - p 2 

where Odlogf?!" 1 ) is independent of n. 
(ii) There are constants C\ and C2 such that 

Qn < Ci(l -p + 5) n - Nl e~^°°s s \ (3.17) 

for all n>N 1 = Ni{5). 
Proof. If n > m, then we have from (|3.7|) 

ra-1 

logq„ = - ^ log(l + S+Pj) + logq m . 
j=m 

Using the inequality | log(l + 1) — t\ < Ct 2 , we obtain 

n-l 

logg n = - pj - 5(n - m) + logg ro + E 1} 

j=m 

where the error term E\ satisfies 

n— 1 n— 1 

|Si|<Cr 1 ^(5 + ft ) 2 <^2^P J 2 . (3.18) 

j—m j—m 

The last inequality above holds since 6 -C p < • Here and in the rest of this proof, Ej 's denote 
errors to be estimated. We then have from (|3.3[) that 



n ^ 

log— = -(5 + p)(n-m)-2p V — - — - 



^1 

j— m+1 
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Since A j L _ 1 = x^-j is decreasing in j, we have 



^ IJ-T + wA log T 



j—m-\-l 



1 - A 



— m— 1 



£ T 

j=m+l 



A-3 



A -n-l 



m+1 



1 - A~ x 



-dx 



< 



A~ 



1-A- 



So we have 



lGg = ~ {6 + P){n - m)+ \^A l0g { 1 - 



£ 2 



where the new error term E2 satisfies 



One can see from Q3.5P that 

where 

So, we have 



2p 



1 + £3 



log A 
\E 3 \ < CVS. 



log ^ = + p)(n - m) + (1 + E 3 ) log f — 



which in turn implies 



f- A" 



wherc 



£4 := pm — 5(n ~ m) + E% + E3 log 



I- A 



-m—1 



l-A- 



Notc that 



Pn - Pn+1 = (Pn+1 ~ p) J^-l-l 



so that 



g m e-P" 1-A 



-m—1 



^4 



Pn - Pn+1 Pn+1 ~P A - I 

Since p n /p = (A n+1 + l)/(A' i+1 - f), we have 

' Pn+P 



(n+ l)logA = log 



Pn - P 



and, since log A = 2p + 0(5), it follows that 



where 



1 , f Pn+P , 
pn = - log +E 5 , 

2 \Pn-p, 



\E 5 \ < C{nS + Vs). 



E 2 



(3.19) 

(3.20) 
(3.21) 



(3.22) 
(3.23) 



(3.24) 



(3.25) 
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We then obtain from (|3.24[) 



Q m ~. : e 

Pn-Pn+1 V Pn+PPn+l-P A- I 



1 Pn-P I -A'" 1 ' 1 p F , . 

-Qm - A : e Ei ~ E *. (3.26) 



\]P 2 U ~ P 2 ~ P A ~ l 

Suppose now that m = Nq — 1 and m < n < N\. Then we have E$ = O ( | log <5 1 ~ 1 ) . We will 
show that 

= l + 0(\logS\- 1 ), (3.27) 



Pn+l - P 

1 — /l~ m ~ 1 

gm A l = 1 + 0(1 log ^r 1 ), (3.28) 

£ 4 = 0(| log (T 1 ). (3-29) 



Once we have these estimates, then (i) follows from p.26[) . 
To prove (|3.27p . we first observe that 

Pn-P A n + 2 -l = 



p n+ i - p A n+1 - 1 V A + A 2 + ■ ■ ■ + A"H 
Since A > 1, n > | log<J| and A = 1 + O(VS), we have 



Pn-P 



= (1 + 0(V*))(1 + 0(| log^- 1 )) = 1 + 0(| log<5|- x ). 



To prove (|3.28[) . we use inequalities 

(m + l)s - -m(m + l)s 2 < 1 - (1 - s) m+1 < (m + l)s 

which hold for all s e [0, 1]. Since A^ 1 = 1 - 2p + 0(5), we have 

(m + l)(2p+0(S)) - ^m(m + l)(2p + 0(S)) 2 < 1 - A'" 1 ' 1 < (m + l)(2p + 0{6)). 

Since m = 0(\ log<5|) and p = 0(-j5), we have 

1 - A-" 1 - 1 = 2(m + l)p + 0(S\ log 8\ 2 ). (3.30) 

Note that 

1 1 1 

A~l ~ 2p+0(S) ~2p~ + 

Since q m = m 1 +1 + 0(^J~8) by Lemma l3.1l and m = Nq — 1, we infer that 
1 - A~ m_1 



A- 1 



So, ([3381 is proved. 



+ 0(\/*)) (2(m + l)p + C>(J| log 6\ 2 )) (J- + 0(1) 

l + oGiogSi" 1 ). 
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To prove (|3.29|) , we first estimate E 2 . We have from (|3.9[) that 

n N N± 



E pHcY,$+ E p. 

j=m+l j=N j=N+l 

N 1 Ni 

^Ei + E^ 



j=N ■ 



N 



< c 



^+^1) =0(|log«5|- 1 ). 



On the other hand, it follows from (|3.30[) that 

pA-" 1 - 1 _ p + Q(S\logd\) 
1 - A-'"- 1 ~~ 2(m + l)p + 0{6\lo£ 

So we infer from (|3 . 1 9[) that 



E 2 ^0(\logS\- 1 ). 
Since p = 0(V8), we obtain from (|3.30[) that 

1 - A-™- 1 



1 > 



l-A- 



T >l-A-™- 1 >cV$|log<f|. 



We then infer from ()3,21[) that 



S 3 log 



l-A- 



1 - A-™- 1 

Thus we have from ([3~23]) . (j3~BT1) and ([3321) that 
|D 4 | < |pm|+5(n-m) + |D 2 | 



< cVS|log$|. 



S3 log 



1-A- 



l-A- 



< c (ys\ ]ogs\ + 1 log^r 1 ) < c\ \ogs\-\ 



so (|3.29[) is proved. 

We have from (|3~22l) that 



= 9m e 



1-4- 



So, it follows from (|3T29)l that 



1_ 



c 2 



qm < Cie 

for some constants Ci and C2. Now, (ii) follows from p,12[) . This completes the proof. 



(3.31) 



(3.32) 



□ 



4 Asymptotic behavior of the singular function 



Let 



Rs := { x e R 3 \ (D 1 U D a ) | p(x) < \ \ogd\- 2 } , 
where p is given by (|2.1j) . Note that Rs is a narrow region in between Di and D 2 . Let 



v(x) :=Utt ^ q n )h(x) = ^q n 



|X-Pr 



X + p„ 



(4.1) 



(4.2) 



In this section we investigate the asymptotic behavior of Vw(x) in the region Rs- We obtain the 
following proposition. 
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Proposition 4.1 For x = (x, y, z) £ Rs, we have 

Vt'(x) = 2§+ 2 p{x)2 ((1,0,0) + 0(| log-T 1 )) ■ (4.3) 

It turns out that |c>j,i>(x)| and |<9 z v(x)| can be estimated without much difficulty. In fact, we 
obtain the following lemma whose proof is given in Subsection 14.11 

Lemma 4.2 For x = (x, y, z) £ Rs, we have 

\d y v(x)\ + |0,«(x)| < (x) (l + Iog(l + ^p)) (4-4) 

for some constant C independent of 6. 

Estimates of d x v(x), especially those terms for No < n < N\, are quite involved. Based on 
Lemma 13.31 (i) we compare v which is given as an infinite scries with the integral defined by 

Wo(x) := I (ix-i o,o)| - |x + m,o)|) 7P=f *• (45) 

where (p, 0, 0) be the fixed point of the combined reflections R 2 Ri- We obtain the following lemmas 
whose proofs are given in Subsection 14.21 and 14.31 respectively. 

Lemma 4.3 For x = (x, y, z) £ Rs, we have 

gx"o(x)= 2(5+ 2 p(x)2 (l + OdlogC 1 )). (4-6) 

Lemma 4.4 For x = (x, y, z) £ Rs, we have 

d x v{x) = a^o(x) (1 + 0(| log^- 1 )) . (4.7) 
Proposition 14. II is an immediate consequence of above lemmas. 

4.1 Proof of Lemma 14.21 



We first observe that if x = (x, y, z) £ Rs, then \x\ < 1 + 5 — y/l — y 2 — z 2 and p < | log S\ 2 , and 
hence 

\x\ < <5 + p(x) 2 . (4.8) 
Using notation p = p(x), v can expressed as 

v(x) = Y q„ I 1 = - 1 = 1 . (4.9) 

h> \^{x- Pn ) 2 +p 2 ^(x+p n ) 2 +p 2 ) 

So, it suffices to estimate \d p v\. Note that 

1 1 \ p P 







P W(*-Pn) 2 +P 2 V(x+Pn) 2 +P 2 J _ l(x-Pn) 2 +P 2 ] 3/2 [(* + Pn) 2 + p 2 ?' 

J-x [(t-Pn) 2 + P 2 f' 2 

Therefore, we have 



2 



y/(x~p n ) 2 + p 2 y/{x+ p n ) 2 + p 2 



- 3(3 L [{t-p n ) 2 +P 2 ? dt - 
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By (HU) we have 

(t- Pn ) 2 +p 2 >C(p 2 +pl) 
for some constant C. It then follows that 



(4.10) 



d 



1 



\J{x-p n ) 2 + p 2 \/{x+p n ) 2 + p 2 



<c p ^ + 5 )<c p 



p 4 + pi 



p 2 + p 2 n 



So we have 



n=0 ' ^ Pn 



Let N = N(S). Using Lemma I3TTI we have 



n=0 



P 2 + 



■f (1/n 2 + p 2 ) n 



< 



pn 



p 2 n 2 



<C 1 



/OS 



1 + p 2 s 2 



< C- 



1 



1 + y 



/9 + Vs 

If n > TV, then q„ = (^(a/^), and thus we have from (13.121) 



C 



^ p 2 +pI ~ ^irP 2 + 6 ~ p + Vs 



This completes the proof. 

4.2 Proof of Lemma [4.31 

Let 



Aiogsr 1 




d x v (x) = / + 


1 d * ( 




J|log5|-i \ 


~I + 11. 





dt 



□ 



If | log < t < 1, then |x± («, 0, 0)| > Ct for some constant C and for all x g R 3 \ (£>j U L> 2 )- 

< C| log5| 2 . 



Since p = O(VS), we also have \Jt 2 — p 2 > Ct. Thus we have 

1 



\II\ < C 



\lo g 5\ 



-i t 3 



(4.11) 



Suppose now that p < t < | log <5 1 1 . Using (|4.8|) and the fact that p = 0(y/~5) again, we have 
for all x £ R$ 



\tx\<t(p 2 + S)<-^-(t 2 + p 2 ), 
\\ogS\ 

M 2 < (p 2 + 5) 2 < (* 2 + p 2 ) 



(4.12) 
(4.13) 
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for some constant C independent of 5. Thus, we have 

1 1 



x±(i,0,0)| 3 (( x ±t)2 +p 2 f /2 {t 2 +p 2 ±2xt + x 2 ) 3/2 

'l + odiog*!- 1 )). 



(f + p 2 ) 3 / 2 

From the mean value property, we have 
-1 1 



|x-(t,0,0)|3 |x+(i,0,0)| 3 
1 



< 



((t 2 + x 2 + p 2 ) - 2xt) 3/2 ((t 2 + x 2 + p 2 ) + 2xtf /2 
6\xt\ 



' \{t 2 + x 2 + p 2 ) - \2xt\f 2 ' 
It then follows from (|4TT2"|) and (grig} that 

-1 1 



|x-(i,0,0)| 3 |x+(i,0,0)| 3 



< CllogcJI" 1 



(t 2 + p(x) 2 ) 3 / 2 



for some constant C independent of S. 
Since 



0.r 



1 



|x-(t, 0,0)| |x+(t, 0,0)| 

x — t X + t 



t 



x-(t,0,0)| 3 |x+(t,0,0)| 3 
1 1 



|x-(i,0,0)| 3 |x+(t,0,0)| 3 
we obtain from (|4.14|) and (|4.16j) 

1 1 



-1 



3e 



|x-(i, 0,0)| |x+(t,0,0)|y [t 2 +p 2 f/ 2 



|x-(i,0,0)| 3 |x+(i,0,0)| 3 



(2 + OOlogcT 1 )) 



It then follows that 



/= (2 + o(i log ^r 1 )) f 

J V 



(t 2 + p 2 f 2 V^—P 1 



(It. 



Using the substitution t = ^/t 2 — p 2 , one can see that 



l log5 l _1 t 1 (\\ogd\- 2 -p 2 ^ 



(t 2 + P 2 f 2 jW—f P 2 + P 2 \\\og6\- 2 + p< 
Since p = y/2S + 0(8) and p < | log<5| -2 , we have 



|log5|- 2 -p 2 ^ 1/2 



and hence 



p 2 + p 2 \ | log 5 1 2 + p 2 
2 



1 



26 + f 



25 + p 2 
(1 + OOlogC 1 )), 



(l + Odlog*)- 1 )), 
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which together with (|4.11|) yields 



9 x v (k) = (l + 0(| \0gS\- 1 )) + 0(| \ogS\ 2 ). 



Since p < | log <5| 2 , the above formula can be written as 



This completes the proof. 



□ 



4.3 Proof of Lemma 14.41 

Let A*o = [| log <5 1 ] and N\ = [ g , ^ gS i ) as before. Let 



JV -1 JVi-l 00 , 

«(*)=£ + E + E d * \r^—^-r-r- 



and 



n=0 


ti—Nq n— 


:=5x(x) 


+ 5 2 (x) + 






Qc«d(x) = / + / + 








:=Ji + J 2 +/ 3 - 





|x-(f, 0,0)| \x+(t,0,0)\J 



We first estimate Si, S3, I-y, and 13. There is a constant O > independent of n such that 
|x ± p n | > Cp n for all x. G Rs- So we have from (|3.9[) that 



JVo-l 



\Si(x)\ < E 



n=0 
iV -l 



1 



1 



|x-p„| |x + p„ 



< C V n 2 - < C| log J| : 



n=l 



We also have from Lemma 13.31 (ii) that 



is 3 (x)i < J2 

n—Nx 
00 



|x-p„| |x + p„| 



< 



00 ^ 

E ^ 



n=N± 



- C E 7 (l-P + <*)" _JVl e Ga ^l*»««l < C. 



n=N! 



Similarly, we have 



|/i(x)| < 



|x-(t, 0,0)| |x+(t,0,0)|; 
<C [ ^dt<C\\og5\ 2 , 



dt 



15 



and by Lemma 13.21 (iv) 
\h(*)\ < 

< C 

So far, we showed that 
We set 



PN 1 



1 



1 



|x-(i, 0,0)| |x+(t,0,0)|; JW^p 



dt 



PN 1 J 



-.dt < C- 



1 



5 2 (x) 



<S 5 /4^tTT^ - 5 S/4 e l/(2V5|log«|) 

|5i| + |5 3 | + |/i| + |/ 3 |<C|log5| 2 . 

1 1 



< c. 



(4.18) 



Ni-X 
n=JV 



|X- P „| |X + P„|7 ^2 _ p 2 ' 



and shall prove 



S 2 (x)-7 2 (x) 



< C 



1 



(4.19) 



Let us first show that Lemma POl follows from (|4.18[) and f|4. 19[) . We observe from Lemma 15751 
(i) that 

5 2 (x) = 5 2 (x) (1 + 0(| log SI' 1 )). 



So we have 



d x v = Si + S 2 + 5*3 

= 5 2 (l + 0(|log ( 5r 1 )) +5i + 5 3 

= I 2 (1 + 0(| log^r 1 )) + (5 a - / a ) (l + 0(| log*)" 1 )) + S x + S 3 
= d x v (l + O(\logS\- 1 ))+R, 



where 



R = -{h + h) (1 + 0(| log^r 1 )) + (5 2 - h) (1 + 0(| logdT 1 )) + Si + 5 3 . 
Since p < | log (5| " 2 , one can see from (|4.18[) and (|4. 19[) that 

1 

+ 



l-RI < C |lo g 5| 



Thanks to (|4.6[) . we now have 

= (i + o(| log^- 1 )) + r = d x v (i + o(| log^r 1 )) , 

which we aim to prove. 

The rest of this subsection is devoted to the proof of (|4.19[) . For No < n < N\, let 



7n(x) 



d x 



|x-(p„, 0,0)| |x+O„,0,0)|; v^TT^ 



Pn - Pn+1 



0. r 



Pn+1 



\x-(t, o.o)! |x + (t,o,o)i; y^rrr^ 



dt 



Let 



/(*) == d x 



|x-(t, 0,0)| |x+ (4,0,0)1; 
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By the mean- value property there is t n € [p n+ i,p n ] such that 



f(Pn)(Pn-Pn+l)- / f(t)dt= {p n ~ Pn+lY 



So we have 



7n(x) < - 



dtd x ( ; 7- 

\|x- {t 

d x 



1 



(t, 0,0)| |x+(t,0,0)| 
1 1 



t=t„ 



1 



|x-(i„, 0,0)| |x+(i„, 0,0)| 



t n 



( 4 2_p2 )3 /2 



(Pn -P«+l) 2 
(Pn - Pn+l) 2 



1 



(7„i(x) + 7ll2 (x)). 



Using (|4.8[) . one can show that 

|x±(t„,0,0)| 2 >C(p 2 +t 2 ), xei?5 
for some constant independent of n. So we have 



7nl 



< 



1 



and 



7n2 



< 



p 3 + ^ V^Fp 2 " 



(p n -Pn+lf 



(Pn -Pn+l) 2 - 



(4.20) 



(4.21) 



(4.22) 



p?+tl (i 2 -f> 2 ) 3 / 2 ' 

If n < N — [-^tj], then we have t n m 1/n and \p n — p n +i\ < C/n 2 by Lemma 13.11 and 
Pn — P > C/n by Lemma 13.21 (iii). So, we have 



N N 

£ 7m<CV , - L 



1 w 1 

n— 1 



/9 3 7l 3 



Note that if p < then 



while if p > y/5, then 



So, we have 



N 1 

y — - 



N 



1 1 1 C 

> < — > — < — ■ 

' 1 + p 3 n 3 p 3 ^— ' n 3 p 

71—1 n— 1 



JV 



Z - 



c 



n=N 



Vs- 



If N < n < N 1; we have from ([3~3]) that 



< p„ - p„+i = 2p 



A- n - x (i- a- 1 ; 



(l- J 4-»- 1 )(l-i4-'»- 2 ) 
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and by Lemma 13.21 (ii) , p n — p Since p n > p for all n, we have 

iVi-l iVi-l 



lnl<C ^2 1 {Pn-Pn+lf 

n=N+l n=N+l {p 3 + p 3 )\ hl +1 - P 2 



oo 

< r? V* - ft 2 A~ 2n 

fgS/2 + p 3) S l/2 A -n/2 

oo 

<C V - § 2 A -(3/2)n 

n=N+l v ' ' 



Thus, we have 



Similarly one can show that 



^ c 

n=N P+^ S 



E 7 " 2 - ^wr 



n=AT 

This completes the proof of (|4.19[) . □ 



5 Proofs of Theorem 11.11 and Theorem 11.21 

Proof of Theorem [LT1 It is helpful to recall that e = 25. We have from (fL2"0)) and Q2TTT} that 

u\dD 1 - u\ dD2 

u = u in — h + 9 

= ~( u \dD 2 - u\ aDl )(4ir^2q n )h + g 



n=0 



= 2( u \dD 2 - u\ dDl )v + g. 



(5.1) 



We emphasize that \Wg\ is bounded on any bounded subset of IR 3 \ (D\ U D2) regardless of S as 
explained in Introduction. Since h is constant on dD\ and dD^., one can see from p.21[) and (11.221) 
that 



u\dD-, - u\ dDl = - [ H ^~ dcr 

Jd(D 1 UD 2 ) vv 



a(Diur> 2 ) 



ov 



a(niun 2 ) 



dH 
dv 



hda 



1 



Eoo 
n=0 Ii 



J2 



D1UD2 



ff(x)— (r(x-p„)-r(x + Pn )) 

Ov 



dH 



- — (r(x-p„)-r(x + P „)) 



da 



C 



H 



2 E^=0 Ir, 
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It then follows from p,15[) that 



l\8D 2 -u\ 0Dl = A| + log S]- 1 )), 



where 0{ \ log <5| ) is independent of H. So we obtain from (|5.1 

C' 



and from (14.3 



Vw(x) 



C 



log<5|(2J + P (x)2) 

and hence (|1.12j) is proved. This completes the proof. 
Proof of Theorem P By (|3lj]) . we have for n < N = N(S) 



((l,0,0) + O(|log ( 5r 1 )) + V 5 (x), 



g n (fT (p„) - H(-p n )) - 
1 



1 



if 



< 



71 + 1 



ti + 1 \ \n + l 
H(p n ) - H(-p n )\ 



1 .O.oW(- 1 



71+ 1 



0,0 



1 



71 + 1 



H(p n ) - H 



1 



n + 1 



0,0 



H{-p n ) - H - 



71+1 



,0,0 



< [ a/5+ — !— ) + CV^— !— . 

71 + 1 / 71+1 



So we have 

JV-l 



N - 1 N 1 / /-, \ / 1 

£ ?n (fT(p„) - i^-p*)) - £ - # -0,0 -if — ,0,0 



)i=0 



(5.2) 



(5.3) 
□ 



< Cy/S\ \og6\. 



On the other hand, since p n is decreasing, it follows from (|3.15[) that 



J2 q n (H(p n ) - H(-p n )) 



n=N 



<C PN J2ln<cVS\\ogS\. 



n=N 



We also have 



n=N+l 



n=N+l 



Combining above estimates, we obtain (|1~1~6|) . 

The formula (|1.18|) is an immediate consequence of (|1.12[) and (|1.16l) . This completes the proof. 

□ 
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